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THE AERODYNAMIC DESIGN OF WINGS WITH CAMBERED SPAN 
HAYING MINIMUM INDUCED DRAG 

By Clarence D. Cone, Jr. 


SUMMARY 

The basic aerodynamic relations needed for the 
design of wings with cambered span having a mini- 
mum induced drag at specified flight conditions are 
developed for wings of arbitrary spanwise camber. 
Procedures are also developed for determining the 
physical wing form required to obtain the maximum 
value of lift-drag ratio at cruise , when the wing 
spanwise camber-line and section profiles are speci- 
fied, by optimizing the wing chord and twist distri- 
butions with respect to both profile and induced 
drags . 

The application of the design procedure is illus- 
trated by determining the physical wing form for a 
circular-arc spanwise camber line. The efficiency 
of this cambered wing is compared with that of an 
equal-span, flat wing of elliptical planform which 
satisfies the same set of flight operating conditions 
as does the cambered wing. 

The wing pitching- moment equations for opti- 
mally loaded cambe red -span wings at design flight 
conditions are also developed for use in trim analyses 
on complete aircraft designs. 

INTRODUCTION 

The theoretical considerations of reference 1 
indicated that the reduction in induced drag- 
attainable with optimally loaded, non planar lifting- 
systems might allow some gain in wing aerody- 
namic efficiency (lift-drag ratio) over that of flat 
wings when the lateral span of the wing system is 
limited by operational requirements. In particu- 
lar, it was shown that substantial increases in 
effective aspect ratio (as compared with that of 
flat wings of equal projected span) could be 
obtained with relatively simple wing forms having 
cambered spans, that is, wings in which the span 


has curvature in a plane perpendicular to the 
direction of flight. Reference 1 pointed out, 
however, that the practical realization of a net 
efficiency increase would depend in considerable 
measure upon the ability to construct cambered 
wings with profile-drag coefficients and structural 
weights which would not greatly exceed those of 
the equivalent flat-span wing producing the same 
lift. Otherwise, the benefits of the increased 
effective aspect ratio attainable with camber 
might be cancelled by increased profile drag, and 
an increased operational lift force. 

The aforementioned effective-aspect-ratio prop- 
erties of camber ed-span airfoils were derived by a 
theoretical treatment of the ideal vortex systems 
necessary for minimum induced drag, without 
consideration of the physical wing forms needed 
to obtain the optimum distributions of vorticity. 
In order to investigate the overall performance 
gains obtainable with such airfoils, the details of 
the physical wing design must be considered. 
These details include not only aerodynamic factors 
but also the effects of the wing structural weight 
on t lie operational efficiency (ref. 1). In general, 
however, structural weight effects can be taken 
into account after the aerodynamic design of the 
wing is specified. The purpose of this analysis, 
consequently, is to develop the general relations 
needed for the aerodynamic design and efficiency 
evaluation of airfoils which have arbitrarily speci- 
fied camber forms (such as circular-arc spanwise 
segments and semiellipse curves) and which possess 
minimum induced drag at specified flight condi- 
tions, without regard to structural weight effects. 
The methods by which the optimum lift loading, 
corresponding to the minimum induced drag of a 
given spanwise camber line, can be determined 
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are presented in reference 1. The relations are 
developed in a form which permits a direct 
efficiency comparison to be made between an 
elliptically loaded flat span wing and an optimally 
loaded (minimum induced drag) wing of arbitrary 
camber having an equal projected span and pro- 
ducing the same lift force. Thus, in addition to 
prescribing the wing shape necessary for minimum 
induced drag with a given ■ camber form under 
specified flight conditions, the results also show 
the relative efficiency of the complete cambered 
wing. The results may be used to investigate 
the maximum values of lift-drag ratio attainable 
with any given camber form for any prescribed 
set of flight operating conditions. 

Although the intention of this paper is to discuss 
only the general design procedure for wings with 
arbitrary camber, without consideration of the 
comparative efficiencies of any specific camber 
shapes, actual calculations are presented for a 
wing having a circular-arc camber line, in order 
to illustrate the design procedures and to indicate 
the order of magnitude of the lift-drag ratio 
attainable for a particular set of operating 
conditions. 

Consideration is given herein only to the design 
of the basic cambered wing. In the design of a 
complete aircraft, the entire vehicle configuration 
must of course be considered and the parasite 
drag of the fuselage and other components intro- 
duced into the design calculations. However, the 
procedures outlined herein specifically for wing 
design are easily adapted to a complete aircraft 
design by proper inclusion of all the component 
drag coefficients. In addition, the relations dis- 
cussed apply equally well in principle to the 
design of more complex wing forms such as those 
having compound camber, or branched tips (ref. 
1), as well as to the design of low-drag hydrofoils. 
Because of the restricted allowable span length of 
most hydrofoil systems, the use of cambered-span 
hydrofoils appears advantageous as a means for 
reducing both the induced drag and support-strut 
drag of conventional systems. The curvature of 
such hydrofoils might also allow some reduction 
in free-surface effects. The theoretical basis and 
practical application aspects of cambered-span 
wings are discussed in reference 1. 

When the use of cambered-span wings is con- 
sidered for specific aircraft designs, a knowledge 


of the pitching-moment characteristics of such 
wings is of considerable importance. The curva- 
ture of the span results in a vertical distribution 
of the drag forces and this distribution in turn 
affects the wing pitching moment. A quantita- 
tive knowledge of this effect is necessary if the 
aircraft is to be designed to possess desirable 
longitudinal trim characteristics. Consequently, 
the subsequent section “Determination of the 
Wing Pitching Moment” presents the develop- 
ment of the various aerodynamic relations needed 
for determining the pitching moment of optimally 
loaded cambered-span wings for any arbitrary 
pitch-axis location. These relations are essen- 
tially confined to the calculation of the pitching 
moment of the optimum wing form at specified 
cruise conditions. 
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SYMBOLS 

aspect ratio 



wing span of flat reference wing 
projected wing span of cambered- 
span wing 

wing total-drag coefficient 
wing induced-drag coefficient 
wing profile-drag coefficient 
section drag coefficient 
wing lift coefficient 
wing lift coefficient for ( LID) max 
section aerodynamic-force (lift) co- 
efficient 

section aerodynamic-force coeffi- 
cient for wing ( L/D) mar 
wing pitching-moment coefficient 
about arbitrary axis through 
point P 
wing chord 

mean aerodynamic chord of flat 
wing- 

wing root chord 
wing total drag force 
wing total drag force intensity 
(force per unit length) 
wing induced drag force intensity, 
pwT 

span wise-camber depth (fig. 3) 
aerodynamic-force loading intensit} r 
(fig- 7) 
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™ (zj 

J 1 r 

- - sec rely 

-1 A o 

flight altitude 

induced drag efficiency factor 
wing lift force 

lift-force intensity (force per unit 
length) 

wing pitching moment 
section pitching moment 

constant, -^—r 

c ° B 2 

, , To/w 0 

constant, 
dynamic pressure 

radius of curvature of circular arc 
wing area of flat reference wing- 
wing area of cambered-span wing, 
r n 

c{$)ds 

J -s t 

arc-span coordinate (fig. 2) 
arc-tip coordinate 
flight velocity 
weight force 

effective down wash velocity 


W 0 

down wash velocity at center of 
span 

X 

Cartesian coordinate in free-stream 
direction 

Xp 

longitudinal coordinate of arbitrary 
axis 

y , * 

Cartesian coordinates of span wise 
camber line z(y ) 

y f 

dummy variable of integration 

z' 

distance of line-of-action of D Q r 
above I r -axis 

Zp 

^-coordinate of arbitrary axis 

Z T 

coordinate of trim axis 

a 

geometrical angle-of-attack function 

a' 

effective angle of attack 

cP* 

effective angle of attack for 
(Z// D) max 

j 

OLi 

induced angle-of-attack function 


geometrical angle of attack of wing- 
root section 

/3 

camber factor, 

7 672 


r 

r 0 

7 

5 

e 

p 

Psl 

(7 

T 

* 


wing circulation distribution 
circulation of wing root section 


nondimensional coordinate, 


V 

b '!2 


nondimensional coordinate, 


2 

672 


wing twist function 
atmospheric density 
atmospheric density at sea level 
atmospheric density ratio, p/p si 

slope of camber-line tangent, tan -1 ^ 
span ratio, 6/6' ^ 


Subscripts: 

C 

condition at cruise 

c/4 

about quarter-chord point of section 

/ 

flat-span wing- 

L 

condition at landing 

max 

maximum 

0 

in plane of symmetry 

P 

arbitrary axis location 

r 

trim-axis location 

Y 

moment about F-axis 

Parentheses ( ) are used throughout to denote 


the independent variable of the function preced- 
ing the parentheses and should not be interpreted 

r 

to mean multiplication. Example: — (*) denotes 

t 0 

r 

that — is a function of ,s*. 

T 0 


THE DRAG POLAR OF CAMBERED-SPAN WINGS 


In reference 1 the relation between the lift 
coefficient and the induced drag coefficient for 
any cambered-span wing is expressed as 


n — 

( ' D ' i ~irkA 


( 1 ) 


where the lift coefficient (\ and aspect ratio A 
are based on the span length and wing area of a 
flat, elliptical-planform wing arbitrarily chosen as 
a basis for determining the induced drag of the 
cambered wing. The factor k is an efficiency 
factor which determines the maximum effective 
aspect ratio kA of the curved wing and is shown 
to depend only upon the geometrical form and 
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length of the camber line when the wing possesses 
the optimum spanwise distribution of circulation 

— ( s ) corresponding to minimum induced drag. 

r o 

The value of k is given by 


k= 


1 r, 

mp 2 b 


M>r jr 

72 J^To 


dy 


( 2 ) 


In this equation \p is the ratio of the span b of the 
reference flat wing to the projected span b ' of the 

curved wing, and 7 =^^ is the nondimensional 

abscissa coordinate of the camber line (fig. 1). 

Thus, if a flat-span, elliptical-planform wing of 
span b and area S is arbitrarily specified, the 
corresponding lift and drag coefficients of any 
cambered-span wing can be determined and the 
total drag polar of the curved wing can be ex- 
pressed as 


C D =C 1}0 - 


irkA 


(3) 


where C Do is also based upon the wing area S. In 
effect, then, the cambered wing is aerodynamically 
equivalent to a flat wing of area S, aspect ratio 
kA, and profile-drag coefficient C Do . The value of 
the profile drag coefficient C Do of course depends 
upon the physical surface area and form of the 
cambered wing at the corresponding lift condition. 
In the following considerations, the factor \p is 
assigned the constant value 1.0, which means 
that the coefficients of the cambered wings are 
being based on a flat wing of equal span, b = b f . 
The significance of the factor p is fully discussed 
in reference 1. 

z 



1 ^ 

Figure 1. — Geometry definitions of flat and cambered- 
span lifting lines. 


In equation (2), k is a dimensionless constant 
for the optimum lift distribution and each of its 
three main factors is also dimensionless, so that 
the expression for k may be simplified to the form 


where 



(*= 1 . 0 ) 



(4) 

(5) 

( 6 ) 


The factors N A and B are constants for any par- 
ticular camber-line shape and apply only to the 
condition of optimum circulation loading for mini- 
mum induced drag. The factor B depends upon 
the form of the optimum circulation distribution, 
and N a depends upon the circulation-downwash 
ratio, which is also purely a function of the arc 
curvature. 

By use of equation (3), the aerodynamic effi- 
/ L 3/2 \ 

ciency parameter L/D l and —jy Jean be established, 

within the limitations of linear airfoil theory, for 
any cambered-span wing. This efficiency index 
can then be directly compared with that obtain- 
able with optimally loaded flat-span wings. 

PROPERTIES OF CAMBERED WINGS HAVING 
MINIMUM INDUCED DRAG 


THE OPTIMUM CIRCULATION DISTRIBUTION 

For any curved wing form having a spanwise 
camber line specified by the symmetrical, dimen- 
sional function z(y) (fig. 2), the span must possess 

r 

a specific circulation distribution - — (s) if the 

F o 

induced drag is to be a minimum at a given lift. 
This optimum circulation distribution can be 


z 



Figure 2. — Cambered-span lifting line. 
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determined by either the conformal-transforma- 
tion or the electrical-analog technique discussed 
in reference 1. These methods can also be used 
to determine the value of N A (eq. (5)), and inte- 
gration of the nondimensional circulation loading 
T/T 0 determines the value of B (eq. (6)). For 
families of camber lines such as circular-arc seg- 
ments or semiellipses, the optimum distribution 
T/T 0 , and the constants N A and B can be expressed 
as functions of the camber factor 13 which describes 
the particular arc members. (See fig. 3.) The 


optimum nondimensional circulation loadings and 
corresponding values of N A and B are presented 
in figures 4 and 5 as functions of /3 for circular and 
semielliptic arcs, respectively, in order to illustrate 
the loading variation with camber for simple arc 
forms. The distribution r/r 0 for the case £=0 
in each figure corresponds to the loading for a flat 
wing, and it is evident that the effect of camber 
is to increase the relative loading of the outer 
portion of the projected span. 

THE EFFECTIVE ASPECT RATIO 


2 



When the values of N A and B have been deter- 
mined for a given camber line z(y ), the value of 
the efficiency factor k can be calculated by equa- 
tion (2), for ^=1.0. Since the effective aspect 
ratio of a cambered-span wing is then equal to 
kA, where A is the aspect ratio of the reference 
elliptical wing of equal span, the induced drag of 
the curved wing will be less than that of the flat 
wing for equal lifts when k^> 1. This result, of 
course, is dependent on the assumption that the 
curved span is optimally loaded at all values of 
the lift coefficient. The variation of the span 
efficiency factor k with degree of camber (3 is 



Figure 4. — -Variation of the optimum nondimensional circulation distribution and efficiency constants of circular-arc 

camber lines with degree of camber. 


686-793 64 —2 
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Figure 5. — Variation of the optimum nondimensional circulation distribution and efficiency constants of semiellipse 

camber lines with camber factor 0. 


presented in figure 6 for the complete families 
(O^S/3^1) of circular-arc segments and semi- 
ellipses for the condition ^=1.0. 

The following considerations of cambered-span 
wings are divided into three parts. The first 
part presents a development of the relations 
necessary for designing wings which will possess a 
maximum value of LjD (for the imposed structural 
restraints) for specified flight conditions. The 
second part presents actual calculations for a wing 
with circular-arc camber as an illustration of the 
design procedure, and as an indication of the 
magnitude of the efficiency obtainable with cam- 
bered wings. The third part presents relations 
for determining the wing pitching moments at 
design flight conditions. 

THE DESIGN OF CAMBERED WINGS 

Consider first a flat reference wing of elliptical 
planform which is specified by stating its span b, 
wing area S, and constituent section profiles. 
The aspect ratio is then also determined, A=b 2 /S. 
The force coefficients are given by 

(7) 

fpWS 


Cd,/— Cdo,/+^| (S) 

These equations in turn can be used to determine 
the wing efficiency parameter LID as a function of 
C z, . 

Consider now an optimally loaded, cambered- 
span wing of projected span length b' equal to that 
of the flat wing and producing an equal lift at the 

same dynamic pressure ~pV 2 . This cambered 

wing can have any desired combination of wing 
twist, chord distribution, and section profile varia- 
tion, subject only to the requirement that at the 
prescribed flight conditions under consideration 
the wing possesses the optimum circulation load- 
ing r(s) for which the maximum efficiency factor 
k applies. If the wing area S of the reference flat 
wing is used as the basis for determining the force 
coefficients of this cambered-span wing, there 
results 

(\= r L (9) 

l ^ 

C D =C D0 +^j ( 10 ) 


I 
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Figure G. — Variation of the efficiency factor k of circular and semiellip.se arcs with the camber factor p. 


These coefficients now determine the lift-drag 
ratio LjD for the cambered wing. 

A comparison of equations (8) and (10) indicates 
that the cambered-span wing can attain a higher 
total aerodynamic efficiency than the flat wing at 
equal lift values if &>1.0, provided of course that 
C Do is not proportionately larger than C DOtfm 
Because C Do depends critically upon the wing- 
section profiles and physical wing surface area, 
the ability to realize the full induced-drag efficiency 
gains offered by span wise camber obviously rests 
upon the possibility of constructing curved wings 
with relatively small surface areas or low drag- 
profiles or both. As is discussed subsequently, a 
considerable latitude exists for designing wings 
that possess low values of C Do at specified flight 
conditions, since the value of k depends only upon 
the attainment of the optimum circulation loading 
T(s) and is independent of the physical means for 
obtaining this distribution. 

It should be noted that the use in this analysis 
of a given flat wing as the basis for determining 
the coefficients of the cambered wing is purely 
arbitrary and is valuable primarily for efficiency 
comparison purposes. In reality, any area S 
could be used to define the force coefficients of the 


cambered wing. The induced drag of a flat, 
elliptically loaded wing is of course independent 
of the physical wing area and hence of the physical 
aspect ratio for a given span and lift, as can be seen 
from the fundamental relation 


This fact is discussed in some detail in reference 2, 
pages 32 and 33, and in appendix A of this report. 

DETERMINATION OF THE WING SHAPE FOR MAXIMUM LjD 

The first step in the design procedure is the 
specification of the nondimensional camber func- 

tion 8(y) where 8= — and y = This func- 

o j £ o j Z 

tion then becomes the corresponding dimensional 
camber line z(y) when the projected span length 
b ' is given. Obviously, the particular value of b' 
to be used for a given design depends entirely upon 
criteria associated with the requirements of the 
specific aircraft mission, and hence no general 
procedure can be stated for determining its selec- 
tion. Once the camber line 5(y) or z(y) is specified, 

r 

the values of k , N A , and B y and the function — (s) 

to 
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corresponding to minimum induced drag can be 
determined. 

The next step is the selection of the area S of the 
flat reference wing of span b ( = 6') to be used as a 
basis for the coefficients. The choice of equal span 
lengths as a basis for comparison rests upon the 
fundamental importance of the span in determin- 
ing the induced drag. (See eq. (Al) and ref . 1.) 
While the value of the area S is purely arbitrary, it 
is convenient to use the value of S corresponding 
to that of the "optimum” flat wing which best 
satisfies the requirements of the specific mission 
for which the curved wing is being designed. Then 
the aerodynamic efficiency of the resulting curved 
wing can be directly compared with that of the 
optimum flat wing. 

Specification of design requirements. — In order 
to design the wing for practical flight operations, 
certain basic design requirements must be speci- 
fied. These requirements will determine, to some 
extent, the design procedure. For the purposes of 
this analysis, however, the design requirements are 
taken as the speed V L corresponding to the desired 
landing speed, the load W L corresponding to the 
total weight supported by the wing at landing, the 
design cruise speed V c , and the initial cruise weight 
TF C . These criteria are of course arbitrary but 
serve as a basis for development of the general 
design relations. Since the landing speed in gen- 
eral is low compared with the cruise speed (depend- 
ing upon the particular aircraft mission), the land- 
ing speed then determines the minimum physical 
area of a cambered-span wing. Under certain 
conditions, take-off speed and weight rather than 
landing speed and weight may be the governing 
factors. 

Determination of minimum value of chord 
function. — A cambered-span wing is usually 
designed to possess its maximum value of L/D at 
some specific cruise flight condition, and the physi- 
cal form of the wing is therefore such that the 
optimum circulation distribution for minimum 
induced drag is obtained. That is, the wing is 
designed primarily for optimum efficiency at cruise. 
It is assumed for this analysis, however, that by a 
proper use of wing flaps or other variable chord- 
wise camber devices, the optimum nondimensional 

r 

circulation distribution - — (,s*) can also be obtained 

r o 

(approximately) at landing conditions. Then the 
values of N A and B apply to both basic flight con- 


ditions. Thus the optimum distribution T/T 0 can 
be used in conjunction with the landing require- 
ments to determine the minimum allowable chord 
distribution of the wing. 

In terms of the nondimensional coordinate 
y=ylQ> '/%), the lift force at any flight condition 

L=pV jf- i r(7)dr (12) 

becomes, when equation (6) is used, 

L=pVT 0 ~B (13) 

for the design flight and landing conditions. For 
landing, L=W l and V—V L , so that 

Wr 

r 0 ,L= % (14) 

plV^B 

Here T 0tL is the circulation in the plane of sym- 

r 

metry of the wing at landing. Since — ( s ) can be 

f 1 0 

established by the methods of reference 1, the 
dimensional circulation distribution T(s) along the 
arc span can be determined as 

r t =^- r 0 , L (15) 

I o 

The aerodynamic-force loading intensity F f 
along the camber line (fig. 7) is given by 

F'(8) = pVT(8) (16) 

The corresponding section “lift” coefficient varia- 
tion Ci (actually c t determines the aerodynamic- 



Figure 7. — Aerodynamic-force loading intensity and its 
relation to lift intensity. 
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force intensity) is given by 


F'(s)_ P V r(s) 
iC \pV 2 c{s) 


(17) 


and equations (14), (15), and (17) give the chord 
distribution for the landing condition as 


c(s) 


W L 

1 T7! 7, 

2 PlV L C l,L ~2 B 



(18) 


This equation, upon specification of the function 
Ci >L , determines the necessary chord distribution 
and minimum allowable physical surface area of 
the wing as set by the landing requirements. 
Obviously, the minimum allowable value of the 
“effective” wing area S' is determined by the 
maximum allowable value of c XtL , where 

S"=J_ 8 ' c(s)ds (19) 


It also follows from equation (18) that if each 
section profile is operated at the same value of 
Ci tL , the chord distribution c(s ) is uniquely deter- 
mined by the camber-line geometry z(y) since 

r 

— ( s ) is a function only of the arc shape. 

1 0 

Although equation (18) determines the mini- 
mum value of the allowable chord length which 
still satisfies the landing requirements, this distri- 
bution is not necessarily the optimum one corre- 
sponding to the maximum value of LjD for the 
wing at a given lift coefficient C L , since the profile 
drag is critically dependent upon the actual 
section chord length. Thus, if the wing is to 
operate at the maximum possible value of L/D at 
a given lift coefficient, the chord distribution must 
be optimized so as to have minimum profile drag. 
This optimization procedure is considered next. 
If the optimum chord distribution is less than that 
required for landing, the limiting distribution as 
given by equation (18) must of course be used. 

The general form of equation (18) for any flight 
condition for which the induced drag is to be 
minimum and in which W , p, and V are specified 
is 


w r 

c(s) "i b~ F w 
|pF ! c,^B r ° 


( 20 ) 


Determination of wing profile-drag coefficient 
and optimum chord function. — When a wing 
section profile is selected, the section drag polar 
c d (ci) and the section aerodynamic-force curve 
Ci(a') are specified. The wing profile-drag coeffi- 
cient is then given by 

C Do =^ c d c(s)ds (21) 

where c d is the section drag coefficient correspond- 
ing to the operating value of c x of the section. 
The drag coefficient as given by equation (21) 
does not, of course, include any of the effects 
such as tip-drag and profile-drag changes which 
may occur with the curved, three-dimensional 
wing. If it is anticipated that such effects will 
be significant, an additive term can be applied in 
equation (21) to account for them. Assume now 
that the wing is constructed entirely of similar 
profiles all of which are operating at the same value 
of Ci) the total wing profile-drag coefficient is then 

J*‘ c(s)ds 

=c d ~ (22) 

The total wing-drag coefficient then becomes 

c °=^ f +33 < 23 > 

In this equation the profile-drag term appears as 
a function of c h where c x is the section “lift” 
coefficient corresponding to any given flight condi- 
tion. Since, however, the wing is being designed 
for maximum L/D at cruise, c d must now be 
expressed as a function of C L (the total lift coeffi- 
cient) so that the function ~ (C^) can be deter- 
mined. 

The relation between c x and C L is established 
by substituting the general weight expression 

W=C L \ P V*S (24) 

into equation (20), and after rearrangement the 
desired relation is obtained as 
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Since the factor in brackets is a constant, the 
relation between the section and total wing-lift 
coefficients can be expressed simply as 


where 


Ci=mC L 


s r 
— FY 


s 



(26) 


(27) 


The individual factors T/T 0 and c are of course 
functions of the arc-span coordinate $ but their 
ratio is a constant when each wing section is 
operated at the same value of c x . (See eq. (20).) 
The total wing-drag coefficient can be expressed as 

C D =c d (mC L ) YSlk (28) 

In order to determine the chord distribution 

that will make C D a minimum for a given value of 
C L j it is necessary^ to find the optimum value of 
the root chord c 0 . The optimum chord-length 
distribution at that C L will then be given by 

c(s)=c~- (s) (29) 

I 0 

In order to find this optimum value of c 0 , the wing 
proffie-drag coefficient 

C D0 =c d (mC L ) ^ (30) 


must be expressed as a function of c 0 . 
equation (27), 

S 


From 


(31) 


and substitution of c(s) from equation (29) yields 



=c„ f S ‘ £ (s)ds (32) 

Since the second integral is a constant for any 
particular camber line when it is optimally loaded, 


equation (32) can be written as 


where 




JT 

r. 


Sec rdy 


(33) 

(34) 


Substituting these values into equation (28) yields 
as the drag polar 





(35) 


where it should be noted that the factor in paren- 
theses is the value of c t for which c d is to be eval- 
uated; it is a functional relation and not a multi- 
plication factor. The equation for LjD as a 
function of c 0 , for a constant value of C Li is given 

by 


L 

D 



n 2 

b— 
' ttJcA 


(36) 


Now if the section drag polar of the selected profile 
is used, a plot can be made at constant C L of the 
variation of LjD with c 0 , and the value of c 0 for 
which L/D becomes a maximum determines the 
optimum chord distribution for that operational 
value of C L . When this procedure is carried out 
for a range of C L values, a plot of LjD against C L 
will determine the optimum cruise lift coefficient 
C L * corresponding to the absolute maximum value 
of wing L/D, denoted by (L/D) maX j and the cor- 
responding value of c 0 . Then, from equation (30) 
the profile-drag coefficient at cruise is given b 4 y 

C D0 =c a (mC L *) ~ (37) 


where S' of course depends upon the correspond- 
ing optimum value of c 0 . 

Since the foregoing procedure determines the 
optimum chord distribution not only for the opti- 
mum cruise condition (f^* and (L/D) ma: r ) but for 
the entire lift-coefficient range as well, the wing- 
dimensions for maximum LjD at any other value 
of C h are also specified. 
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It should perhaps be emphasized that the rela- 
tion between LID and C L as given by equation (36) 
for the optimum value of c 0 at each C L is not the 
variation of a particular wing LjD as the angle of 
attack is increased. It merely gives the relation 
between LjD and C L for the optimum wing form 
at a given C L value, and for each value of C L the 
wing form will be different (different chord 
distribution) . 

The optimum cruise altitude. — If cruise at the 
most efficient lift coefficient C L is to be attained, 
the wing must fly at the altitude where the opti- 
mum equilibrium condition 

W c =C L *hj cPsl V ( ?S (38) 


is satisfied, since the cruise speed V c is specified. 
(Alternately, of course, the cruise altitude could 
be specified and the cruise speed calculated. If 
both the altitude and speed at cruise are specified, 
the necessary C L for cruise is determined and the 
foregoing procedure gives the maximum value of 
LjD at this C L value.) The subscript C denotes 
the cruise values of the variables and <j c is the 
density ratio p c /p s i • The density ratio is of course 
a function of altitude h and <j c — <T(hc). The opti- 
mum cruise altitude is then given implicitly by 


(T C = 


w c 


iPs 


,Vc 2 S 


(39) 


through the altitude-density function a(h). The 
air density at cruise altitude is 

Pc = cr c p s i (40) 


effective down wash velocity w(s) is shown in 
reference 1 to be given by 


w(s)=w 0 cos t(s ) (42) 

In this equation w 0 is the maximum down wash 
(occurring at the center of the span) and r is the 
slope of the camber-line tangent at point $. The 
effective down wash at a point s of the camber line 
is the component of the total induced velocity at 
point s normal to the camber line (fig. 8). In the 
case of symmetrical arc forms the total induced 
velocity is w 0 and is constant along the arc. 

From equation (5) 


w 0 



( 43 ) 


and if equations (13), (42), and (43) are used 
with the condition W=L, it follows that 

W 

w(s) ^ rh ,-12 cos T ( s ) (44) 

p v\^\ b n a 


Then from the expression for the camber line, 
z(y)j there follows 

tan t=~ (45) 


whence 



(46) 


This equation defines the relation cos r as a 
function of the coordinate y, and since 


In equation (39), W c is the initial cruise weight 
supported by the wing. In practical aircraft op- 
eration with conventionally fueled engines, W c 
will decrease as fuel is used up and if the aircraft 
is to fly at C L * continuously, the altitude must 
increase in such degree that 

Wn 

— ^Constant (41) 

The effective downwash distribution at cruise. — 

The effective downwash distribution along the 
cambered span at cruise must now be determined 
so that the geometrical twist necessary for opti- 
mum loading of the wing can be calculated. The 


s= 



1 + 



(47) 



Figure 8. — Diagram of induced velocity and its relation 
to effective downwash. 
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cos r is also determined as a function of 5 (the 
arc-span coordinate). Thus, the effective down- 
wash distribution corresponding to minimum in- 
duced drag becomes known for any flight condi- 
tion. For the cruise condition, use of equation 
(44) yields 

H7 

w(s) = j=r£p cos t(s) (48) 

P°Vc |_y j BN A 

Wing twist function for minimum induced 
drag— The induced angle-of-attack distribution 
oii(s) is defined by 


ai(s)= tan“ J (49) 

where w(s) is the effective downwash distribution 
along the arc span. At cruise conditions, sub- 
stitution of equations (39) and (40) into equation 
(48) gives 


w 

Vc 


( B )- P ** 8 

rb'i 2 

2 [t] m ‘ 


cos t(s) (50) 


The induced angle of attack at cruise is, therefore, 


<*i(s)= tan 1 -s 


C L *S 


*[U 


cos t(s) 


f (51) 


BN a 

This equation can alternately be written as 

]2 ' W (52) 


, so C L *S . , 

<^l(s) 2 | r 1 / 12 COS T\S) 


Trk[b'} 2 


by using the usual linearizing assumptions and 
equation (4). 

The required geometrical angle-of-attack dis- 
tribution a(s) at cruise is given by 

a(s) =«*($) +<*'* (53) 


where a'* is the section angle of attack correspond- 
ing to the optimum section lift coefficient c t ~ 
mC L . Under the foregoing assumption of geo- 
metrically similar profiles for all sections, a'* is 
of course constant. Substitution of equation (52) 
into equation (53) gives the geometrical angle of 
attack (measured as the angle, in the plane normal 
to the camber line, between the direction of V 


and the section chord line) at cruise in terms of 
the camber and design lift coefficient C L * 

a(s)=a'*+2 cos t(s) (54) 

If the geometrical twist function 6(s) is defined as 

6(s) = a(s) — a 0 (55) 

where a 0 is the geometrical angle of attack of the 
root wing section (at s=y= 0), it can be seen that 
the twist for an optimum cambered-span wing 
constructed of similar profiles will be negative, 
since a decreases as r increases, as indicated by 
equation (54). This relationship means that the 
wing must possess geometrical “washout” for 
minimum induced drag at cruise. 

Final wing form, — With the determination of 
the twist function, all the information necessary 
for the aerodynamic specification of the total wing 
form for maximum cruise efficiency is available. 
The camber line z(y) specifies the frame upon 
which the chosen section profiles are to be arranged. 
For simplicity, the camber line can be assumed to 
pass through the quarter-chord point of each 
section or else the location of the profiles relative 
to the camber line can be specified. The optimum 
chord distribution c(s), as determined by equation 
(29) when the optimum root chord length is 
specified, follows directly from the optimization 
procedure for maximum L/D. Finally, equation 
(55) for the twist function specifies the angular 
arrangement of the section profiles along the 
camber line, and a 0 gives the geometrical angle 
of attack of the center wing section. 

ADDITIONAL DESIGN CONSIDERATIONS 

Although the design relations presented have 
included the principal factors involved in optimiz- 
ing a wing form for specified operating conditions, 
several other special factors which are of im- 
portance must be considered in a complete wing 
design. These are briefly discussed in this section. 

Section profile selection. — A small effective 
wing area S' is particularly desirable for cambered 
wings, especially when use is made of the thick, 
laminar-flow profiles, such as characterized by the 
NACA 65 3 -618 wing section. Since the section 
drag coefficient of such profiles is very small and 
is essentially independent of the section lift co- 
efficient over a wide lift range (the drag bucket), 
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■4 


use of such profiles operating at a relatively high 
value of Ci and with a small chord length can 
result in a substantial lowering of the wing' profile 
drag at cruise as compared with that of more 
conventional sections, provided that the full ex- 
tent of laminar flow can be realized and that the 
cruise Reynolds number is not so large that the 
drag bucket has become very narrow. In addi- 
tion, the relatively large thickness ratio of these 
airfoils provides considerable depth for housing 
the wing-spar structure. The attainment of a 
small value for S' } however, depends upon the 
maximum value of Ci that can be attained for 
satisfaction of landing or take-off requirements. 
It thus appears highly desirable when such pro- 
files are used to make full use of such high-lift 
devices as multiple-slotted flaps in order to in- 
crease c^l to a value sufficiently high for low-speed 
operation. 

In order to evaluate various wing-section pro- 
files for optimizing the wing (L/D) max value at 
cruise conditions, the function 


Cl 

Cd 


C d 


C l 

S' Cl 2 
S 'irkA 


(56) 


can be optimized for a series of different profiles, 
and the profile having the largest value for this 
parameter can be determined. The final profile 
selection must of course also be based upon the 
Ci t max value and upon a sufficient thickness ratio 
to meet structural requirements. 

The optimization considerations are based upon 
the restrictions that the wing be constructed 
entirely of similar profiles all operating at the 
same value of c z . It is possible, however, to 
design the cambered-span wing so that it will 
possess the maximum possible value of LjD when 
the restrictions of operation at constant c t and 
similar profiles are removed. When every sec- 
tion of the wing is simultaneously operating at 
its own value of (£/ /D') mox , the total wing will 
then have its maximum value of LjD for a given 
lift. For example, if the cruise conditions W c , 
V c , and h c are specified, the value of (L f ID') max 
can be determined for each section by use of the 
relation 

U__Ci{a')qcC c os r 
D'~D t '+c d W) qc c {b7) 


where 


and 


Di ~p c wT c 



c cos r 


L' = pcV c T c cos r 
=Ciq c C cos r 


(58) 


(59) 


Under the given cruise conditions, equation (59) 
leads to the condition 


Cl c = Constant (60) 

for any specific section and use of this condition 
in equation (57) will ultimately determine the 
optimum chord length, angle of attack a, and 
twist angle for each section, corresponding to the 
value of ( L'/D ') max . The section data used 
should of course correspond approximately to the 
anticipated flight Reynolds number of the section. 

Angle-of-attack variations. — The optimum cam- 
bered-span wing is designed primarily for cruise 
conditions, and the geometrical twist provides the 
angle-of-attack distribution necessary for the 
optimum circulation loading at cruise speed. 
For lower speeds an increase in section angle of 
attack is necessary if a constant lift is to be main- 
tained. Because of the camber of the span, 
however, an increase in the effective angle of 
attack of the center wing sections obtained by 
pitching the wing will have less effect on the out- 
board section angles of attack, due to the effect 
of camber, if the wing is rigid. A possible solu- 
tion to this problem lies in the use of tr ailing-edge 
flaps. By a suitable programing of flap deflec- 
tion, the section lift coefficient can be varied 
while maintaining a constant section geometrical 
angle of attack (corresponding to the cruise 
condition). In this way, the wing could always 
operate at essentially the same flight attitude for 
take-off, landing, and cruise. In addition, it 
would be possible to maintain an approximately 
optimum circulation distribution at all flight 
conditions without need for pitching of the entire 
wing. The use of such flaps is also quite com- 
patible with the need for high values of Ci tmax 
and small values of S' at cruise. 

Take-off and landing requirements. — The design 
procedure outlined above is based upon the speci- 
fied landing weight and speed. Inasmuch as the 
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take-off weight is usually considerate greater than 
the landing weight, for the same Ci r7nax the take-off* 
speed must therefore be greater than the landing 
speed. The take-off speed could be equally well 
used as the criterion for determining the minimum 
chord distribution. However, the desirability of 
maintaining a reasonably small value of S' for 
optimum LjD operation may require some com- 
promise of the minimum landing or take-off speeds, 
unless means for obtaining very high section lift 
coefficients, such as external boundary-layer con- 
trol, are utilized. For special missions where 
cruise efficiency is critical, use of jet-assistecl take- 
off is quite feasible as a means for overcoming any 
take-off problems which might exist. 

Structural considerations. — The preceding de- 
sign outline has been based solely upon aero- 
dynamic considerations without regard to struc- 
tural factors. The structural factor of primary 
importance is the wing structural weight as com- 
pared to a flat-span wing producing equal lift. 
In a complete design analysis, the effects of struc- 
tural weight must of course be considered. The 
effect of the span ratio \p on weight properties, as 
well as the possibility of utilizing lightweight 
aeroelastic wing structures for weight reductions, 
is discussed in reference 1. 

Special missions. — In addition to the possibility 
of using cambered-span wings for improvement of 
the efficiency of conventional aircraft, there are a 
number of special areas where application of such 
wings may be of significant value. Since the 
primary advantage of cambered wings lies in their 
high effective aspect ratios for a given span length, 
operational missions where high induced drag is a 
serious problem are of particular interest. There 
are several such missions. In the operation of 
STOL aircraft, the very low take-off and landing 
speeds require operation at very high C L values 
with very high D x . For such airplanes, operating 
conditions often severely limit the allowable span 
length and hence the geometrical aspect ratio. 
In addition to the very high induced drag en- 
countered by such aircraft, it is predicted in 
reference 2 that at very large values of C L (or for 
small values of A ), the airfoil lift can be appre- 
ciably reduced by the large downwash velocities. 
Use of spanwise camber could help in alleviating 
both problems. In situations where the physical 
span length is limited, camber can be used to 


increase the effective value of the aspect ratio 
over that of flat-span wings. Various other areas 
include efficient subsonic cruise at extremely' high 
altitudes (where the low density requires high C L 
values), high endurance aircraft, and cargo aircraft 
operation with increased payload. Possible sta- 
bility and control applications are also suggested 
in reference 1 . 

ILLUSTRATIVE DESIGN OF A CAMBERED-SPAN 
WING 

In order to illustrate briefly the application of 
the foregoing design procedure, calculations have 
been carried out for a cambered-span wing which 
satisfies an arbitrary set of prescribed operating 
conditions. The conditions selected are typical 
of current operational requirements met by con- ■ 

ventional-winged transport aircraft. No attempt 
has been made in this illustration to determine 
the absolute optimum wing form as regards 
(. L[D ) max , and the wing chord distribution used 
corresponds to the minimum value as set by the 
prescribed landing condition. The wing section 
chosen for the illustration is the NACA 65 3 -618 
profile. However, the use of laminar-flow profiles 
is in no way mandatory with cambered wings, 
although such profiles possess some advantages 
because of them low drag coefficients. It is 
assumed for purposes of illustration that the 
section profile-drag coefficients obtained in wind- 
tunnel tests can be realized in full-scale, practical 
operation. In reality, extremely smooth wing 
surfaces would of course be required for extensive 
areas of laminar flow. The results of this study, 
therefore, may be considered as a realistic, al- 
though not conservative, indication of the general 
level of aerodynamic efficiency of the particular 
cambered wing for the given operating require- 
ments. ^ 

It should perhaps be emphasized that the rela- 
tive efficiencies of the flat and cambered wings of 
this illustration are valid only for the particular ^ 

set of operating conditions imposed. The particu- 
lar camber line used was quite arbitrary and use 
of a more optimum camber form might result in 
significant increases in wing efficiency. In any 
specific application, the most optimum camber 
form can only be found by a series of design 
trials based upon the particular operational 
requirements. 
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BASIC DESIGN CONDITIONS 

The cambered wing arbitrarily chosen for this 
illustration has a circular-arc camber line with 

£=0.8 and ~^=58 feet. Thus, when optimally 

loaded, the wing will have the following values 
for its constants: &=1.32, N A — 2.561, and 
1.619, as shown in figures 4 and 6. The basic 
flat wing chosen for comparison purposes and as a 
basis for coefficients, has the following character- 
istics: elliptical planform, — — 58 feet, A=8, 

z 

$=1,683 square feet, and the wing section is the 
NACA 65 3 -618 throughout. This same profile 
is used for the cambered wing. 

The prescribed speed and weight conditions are 

W £=85,000 lb ^=110,000 lbl 

y (6i) 

^£ = 105 mph Fc = 310 mph J 

These operating requirements and the span and 
wing area of the fiat reference wing are typical of 
the operating specifications of an actual transport 
aircraft wing. 

THE SPANWISE CAMBER-LINE GEOMETRY 

The equation of a circular-arc span wise camber 
line is given by 


the value of r is given by 


1 l+£_ 2 

2 2 £ 


(63) 


and substitution of equation (63) into equation 
(62) gives the dimensional camber line z(y). For 

b' 

the present configuration where £=0.8 and 58 

feet, the function z(y) is as shown in figure 10. 

Since most of the wing design variables are 
expressed as functions of arc length s instead of y 
(fig. 2), the relation s(y) is needed. This function 
is obtained by use of the arc-length expression 

-JT V i+ GfT *'="• (64) 

which, when applied to equation (62), yields 

s(y)=r sin- 1 ^ (65) 


Equation (65), when plotted for the present case, 
gives the curve shown in figure 11. This curve 
allows direct conversion from one system of 
coordinates to the other. 

It is also necessary to determine the function 
cos t(s). The relation for r(y) is given by 


2 =r---[r 2 — -?/2 ]i/2 

= Kl-COSr] (62) 

where r is the radius of curvature of the arc. (See 
fig. 9.) For any values of j3 and semispan 6'/2, 



Figure 9. — Geometry of a circular-arc camber line. 


r(y)=tan ~‘ % 

= tan' 1 V = 

- \**—y 2 

^sin" 1 ^ (66) 

where r is given by equation (63). Use of equa- 
tion (65) allows the determination of r(.s), and 
the function cos t(s) follows directly. Alternately, 

cos r (y) == ^ 1+y2(r 2 _^p 

=V i ~S ( 67 > 

The function cos r(s) for the present camber line 
is plotted in figure 12. 
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Figure 10. — Circular-arc camber line used in the illustrative calculations. 


0—0-8, 58 feet. 
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Figure 11. Relation between are length coordinate s and span coordinate y for a circular-arc camber line of /S=O.S 
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THE CHORD FUNCTION c(s) 

Although the foregoing procedures allow the 
optimum chord distribution to be determined for 
the cambered-span wing, the chord distribution 
to be used for the purposes of this example 
corresponds to that set by the landing require- 
ments listed previously. Thus, the chord dis- 
tribution used is not necessarily optimum. 

The chord function is obtained by substitution 
of the landing conditions of equation (61) into 
equation (18), with p L taken as 0.002378 slug 

p 

per cubic foot. The function ($) is given 

-L o 

r 

indirectly in figure 4 as — (7), and it is assumed 

1 o 

that by use of a suitable slotted-flap arrangement 
the value c XtL = 3.0 is operationally attainable 
with the NACA 65 3 -6 1 8 profile. Values of 
Ci'L as high as 3.2 are obtainable with double 
slotted flaps (see ref. 3), and the chosen value 
of 3.0 is not unusually optimistic, even without 
the use of synthetic means of boundary-layer 
control. The resulting chord distribution c(s) 
is shown in figure 13. 

THE WING DRAG POLAR FOR MINIMUM INDUCED DRAG 

Use of equation (19) yields the effective physical 
area of the wing S' as 1,340 square feet. Equation 
(27) gives the value of m as 1.675, and the relation 
between total and section lift coefficient becomes 

1.675 C L (68) 


Thus, in conjunction with the section drag polar 
for the NACA 65 3 -6 1 8 section (see fig. 14), the 
wing drag function C D (C L ) follows directly from 
equation (28). The “polar” for the present wing 
series is presented in figure 15, where the function 

(C L ) is also plotted. From this latter curve 
the value of C L for maximum LfD is 

Ci*=0.315 (69) 

The polar of figure 14(a) is for a Reynolds 
number of 9X10 6 , while the cruise flight Reynolds 
number (based on the wing root chord) is approxi- 
mately 27X10 6 . Although experimental data 
on the NACA 65 3 -6 1 8 section does not exist for 
Reynolds numbers above 9X10 6 , there is evidence 
that similar profiles operating at the design c x value 
of 0.528 have even lower values of the drag coeffi- 
cient for Reynolds numbers up to 35X10 6 , (See 
ref. 4, fig. 14, p. 18.) With increasing Reynolds 
numbers above 9X10 6 , however, the extent of the 
low-drag bucket decreases slowly and for high 
values of c x (c*>0.8) the drag coefficient (for 
sections near the wing root) will probably be 
slightly greater than shown in figure 14(a). 
Still, since most of the wing will be operating at 
local Reynolds numbers considerably below 
27X10 6 , where the section drag coefficient is 
actually lower than indicated in figure 14(a), 
the total wing drag coefficient variation should 
be very close to that based on figure 14(a). 
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Figure 14.— Section coefficients for the XACA 603-6I8 
profile. 

THE OPTIMUM CRUISE ALTITUDE 

From equations (39) and (69) there follows 



(b) Lift curve. 
Figure 14. — Continued. 


<r c = 0.845 (70) 

whence 

A c — 5,600 ft (71) 


At this altitude, the Mach number for cruising 
flight is 0.416, which is considerably below the 
critical value at the section operational lift 
coefficient (see fig. 14(c)). 

THE EFFECTIVE DOWNWASH AT CRUISE 

From equation (50) and figure 12 the down wash 
distribution w(s) may be determined, and by use 


of either of equations (49) or (52) the induced 
angle-of-attaek distribution may be calculated. 
These two functions are presented in figure 16, 
for the present wing. 

THE TWIST FUNCTION d(s) 

The value of c/* (eq. (53)), corresponding to the 
optimum section force coefficient c t ~ m(' L *= 
0.528, is obtained from the section lift curve shown 
in figure 14(b) and equals 0.4°. Equations (53) 
and (55) can be used to calculate the wing geo- 
metrical angle-of-attack and twist distributions 


u 


w 


0 
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(c) Critical Mach number. 
Figure 14. — Concluded. 



(a) Drag polar. 

Figure 15. — Total drag polar and LjD variation for the 
cambered wing. 

a(s) and d(s). These functions are plotted in 
figure 17 for the example wing. 

THE FINAL WING FORM 

The relations shown graphically in figures 10, 
13, and 17, together with the specification of the 



(b) LjD variation. 

Figure 15. — Concluded. 

NACA 65 3 -618 section, fully determine the 
physical form of the wing for cruise at the maxi- 
mum value of LID (for the particular chord 
function being used). As can be seen from figure 
13, the variation of the chord length along the 
arc span is different from an elliptical distribution, 
especially along the outer 00 percent of the arc 
span. From figure 17 it is seen that the opera- 
tional angles of attack at cruise are very small, 
due primarily to the high amount of camber of the 
basic section. Only a very slight amount of 
washout is required, the total twist amounting to 
only —0.85°. 

The LjD variation of the flat reference wing is 
plotted in figure 15(b), where the maximum value 
is seen to be 37.2 at a lift coefficient of approxi- 
mately 0.35. The degree to which the efficiency 
of the flat wing can be improved is indicated in 
this same figure by the curve labeled “Minimum- 
area flat wing.” This LjD variation is for a 
flat-span wing which possesses a minimum chord 
distribution based upon the same operational 
conditions as the cambered-span wing. This 
minimum-area flat wing is assumed to have a 
proper twist distribution so that minimum 
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Figure 16. — Down wash and induced angle-of-attack distributions for the cambered wing. 


induced drag is obtained. The effective area of 
this wing is only 1,000 square feet. It is evident 
that appreciable decreases in the profile drag are 
theoretically possible when the physical wing- 
area is reduced. 

The LjD variation for the cambered-span wing 
is shown as the top curve in figure 15(b), where it 
is seen that a maximum value of 50.0 occurs at 
Ci,— 0.315. The cambered-span wing (at least for 
the conditions and assumptions of this example) is 
more efficient than either the reference or mini- 
mum-area flat wings for a given lift. Although 
the profile drag is a significant part of the total 
drag at the value of C L for ( LjD) max that occurs in 
this particular example, it is of much less impor- 
tance at the higher values of C L where the induced 
drag predominates. 

All the comparisons cited are based on the use of 
the NACA 65 3 -618 airfoil section and merely indi- 
cate the relative efficiencies of a specific cambered- 
span wing (circular arc, 0=0.8) and a specific 
flat-span wing under a particular set of specified 


operating conditions. Neither of these wings can 
be said to possess maximum efficiency in the sense 
of absolute values obtainable. It is possible that 
by use of other airfoil sections and other spanwise- 
camber forms, and with more optimum chord- 
length distributions, even higher efficiencies could 
be obtained with both the flat and cambered-span 
wings. 

An additional point of interest in figure 15(b) is 
the value of LjD for C L values greater than C L *. 
Even when wings are designed for operation at C L 
values considerably above C L * } the Z/Z) of the 
cambered-span wing is still considerably greater 
than the ( L/D) max of the fiat-span wing, at least for 
the assumed conditions of this specific example. 
This difference means that with the specified value 
of W c , the cambered wing could cruise much more 
efficiently at high altitudes than could the flat 
wing. As C L increases beyond 0.7, however, the 
cambered wing becomes less efficient than the basic 
flat wing. This efficiency loss, of course, is purely 
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Figure 17. — Geometrical angle-of-attack and twist distributions for the cambered wing. 


the consequence of the relatively small section 
chord length of the cambered-span wing, resulting 
from the completely arbitrary choice of a mini- 
mum chord length to satisfy the landing conditions. 
This chord distribution may in no way he the 
optimum distribution for large (\ values. At high 
values of (\ the sectional operational lift coefficient 
Ci = 7nC L exceeds the extent of the laminar flow 
capabilities of the particular section being used, 
and the profile drag becomes exceedingly large (on 
the same order of or larger than the induced drag). 
For flight at high C L values the cambered- wing- 
profile drag can be considerably lowered by in- 
creasing the section chord length to its optimum 
value, by means of the optimization procedure 
previously outlined. This procedure of course 
would also aid in reducing the value of Ci, max re- 
quired for landing, 

DETERMINATION OF THE WING PITCHING 
MOMENT 

The pitching-moment equation for the most 
general case of a wing composed of an arbitrary 


development of section profiles, each operating 
with its optimum value of c t and chord length c is 
considered first. Then the equation for the special 
case of a wing which is constructed of similar 
profiles throughout, with each section operating 
at the same Ci value, is derived. It is assumed 
in both of these cases that corresponding values 
of the functions c t (s) } c d (s), c di (s), c w>c/ 4 (s), c(s), 
and w(s) have been determined by the foregoing 
methods and that the corresponding flight dynamic 
pressure at cruise, y 2 pV 2 , is also known. The value 
of c m ,c / 4 is, of course, determined from the section 
pitching-moment curve c W(C/4 (cO for the c t value 
at which the section is operating. 

The equations are developed for two pitch-axis 
locations. The first location is for the T-axis of 
the wing camber line and the second is for an 
arbitrary axis location. The derivations that 
follow are based upon the assumption that both 
types of wings are so constructed that the camber 
line z(y) passes through the quarter-chord point 
of the section profiles. For wings which are not 
constructed in this manner, proper account must 
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be taken of the displacement distance of the local 
profile quarter-chord point from the spanwise 
camber line. 

PITCHING MOMENT OF THE MOST GENERAL WING FORM 

The section pitching-moment contribution. — 

Consider first the section with the quarter-chord 
point located at the point (; y , z) on the span line 
of a cambered-span wing. This profile lies in a 
plane which is normal to the camber line at (y } z) 
and thus makes an angle r with the vertical, where 
r is the slope of the tangent line at the point. (See 
fig. 7.) At the design operating conditions the 
section will produce a pitching moment of inten- 
sity Me// about the spanwise camber line. This 
moment is a function of c h the operating force 
coefficient of the section, where 

F f 

Ci= T^— (72) 

r V2c 

The component of this moment about an axis 
through the quarter-chord point and parallel to 
the F-axis of the camber line is M c/4 ' cos r. The 
force system creating this moment is shown in 
figure 18 for the particular case where the section 
chord is alined with the direction of the free-stream 
flow. For the coordinate system shown in this 
figure with the origin at the center of the camber 
line, the section pitching moment about the hori- 
zontal axis through the quarter-chord point is 


line-of-action of DJ above the F-axis, being 
positive in the upward direction. The positive 
direction of the resulting pitching moment is as 
shown in figure 18. The pitching moment of 
this section about the F-axis of the airfoil is 
given by 

My = —L'x+Do'z'+D/z (74) 

and since 

D 0 'z'=D 0 '[z'-z+z] (75) 

the moment becomes 

M v '=M cH ' cos t+D'z (76) 

Consider now the case where the pitch axis is 
located at any arbitrary position having coordi- 
nates with respect to the F-axis of the camber 
line x P and z P , as shown in figure 19. The section 
moment about this axis is 

M/=~TAx-~x P ]+D 0 '[z'~Zp]+D/[z-Zp] 

= -L'x+D 0 '[z'-z] + L'xp+D'[z-Zp\ 

=M cn ' cos r+L'xp+D'lz-Zp) (77) 

The subscript P denotes the arbitrary axis posi- 
tion and z P is defined as positive when the pitch 
axis lies above the F-axis of the airfoil. 

When the section force coefficient c t (s) (eq. 
(72)) is specified, along with the chord 
function c($), the pitching moment of each section 
can be determined for either of the axis locations 


M c/ ± cos t=— L'x-\-D 0 f [z f — z] (73) 

where L' = F ' cos r and z f is the distance of the 
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Figure 18. — Force and moment system for a wing section 
at an arbitrary spanwise location. 



Figure 19. — Moment relations for an arbitrary axis 
location. 



THE AERODYNAMIC DESIGN OF WINGS WITH CAMBERED SPAN HAVING MINIMUM INDUCED DRAG 23 


of equations (76) or (77) since 

M c/ 4 '(s)= Cm , c , 4 ( S )^F 2 c ( S ) 2 


line z(y), 


(78) 


1 



( 86 ) 


The section moment, moment coefficient, and 
chord-length distributions are continuous func- 
tions of the arc-length coordinate s. 

Wing pitching moment about F-axis. — The 
total pitching moment of the wing, with regard 
to the F-axis, is obtained by integrating the con- 
tributions of the sections along the arc span, 

My= J_ S ] My'ds (79) 

and substitution of M Y f from equation (76) yields 
M y = f * [M c /i cos t+D'z] ds (80) 


and then as a function of s from the relation 




(87) 


From equation (87), z can also be determined 
as a function of s. Since all factors of the inte- 
grand of equation (85) are variables in the general 
case being considered, this integral cannot be 
further simplified. 

Wing moment about arbitrary axis. — The 
wing pitching moment M P about an arbitrary axis 
P(x P ,z P ) follows directly by substitution of 
equation (77) into the integral 


In this equation the factors are all functions of s. 
The moment M c/ / is given by equation (78), 
where c m>c/4 is a function of the sectional force 
coefficient c t . The drag force D' is made up of 
two components, the profile drag D 0 ' and the 
local induced drag D/ where 


Do —c d qc 

(81) 

and 


Di =F' tan a t 

(82) 

As has been shown previously 


F'=c iq c 

(88) 

and 


ft 

n 

^18 

(84) 


where w is the induced velocity component 
normal to the lifting line (the “effective” down- 
wash, as given by eq. (42)). 

Thus, equation (80) can be written in the 
expanded form 

M Y =q.^ s J^c m _ c/i c 2 cos r + 2 ^c+Ci^cJ J-rfs 

(85) 

The function cos r(s) can be determined as a 
function of y from the geometry of the camber 


M P =J_ S ‘ M/ds (88) 

Thus, 

Mp—qJ 1 ^c mtC/4 c 2 cos.r+^c+c^cj 2 

+XpC ( c cos t— (89) 
Alternately, equation (89) can be written as 
Mp=My+qx P j CiC cos rds 

~1 z rj ‘ (j^+ c * yj c<ls (9°) 

PITCHING MOMENT OF A WING CONSTRUCTED OF 
SIMILAR PROFILES 

The pitching-moment relations for a cambered- 
span wing which is constructed entirely of similar 
profiles, all operating at the same value of the 
section force coefficient c t can now be derived. 

Wing pitching moment about F-axis. — The 
pitching moment about the F-axis of a wing with 
similar airfoil sections along the span can be ob- 
tained directly by specializing the general form of 
equation (85). Under the design restrictions, it 
follows that the chord distribution c(s) can be 
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written as 


c(s)=Co 


n 

To 


(91) 


moment coefficient can be defined on the basis of 
the wing area S and mean aerodynamic chord c 
of a flat- span reference wing 


where c 0 is the wing root chord and V/T 0 is the 
nondimensional circulation distribution for mini- 
mum induced drag. Thus, equation (85) becomes 


C mt p — ; 


M P 

~pV 2 Sc 


(98) 


M y =qc 0 ‘ (jr) cos rds 

+c tf J* ‘ £ zds+c, J 1 ~ y sc/sj (92) 

Pitching moment about an arbitrary axis. — The 

pitching moment about an arbitrary axis P is 
obtained directly by specializing equation (90), 

M P —M Y -\-qc 0 ~CiX P J * cos rds 

r< r , r*t r w , 1 

Since the following relations are valid : 


s t 


r , _b' n 

T a ds 2 G 


W W 0 
7=7 COS r 


(93) 

(94) 

(95) 


f s ‘ r w j w 0 f s < r , 

J-„F.y*“Fj..,r; c08T * 


y w 0 R 
= 2 v B 


(90) 


equation (93) becomes 

Mp=M y -\-qc 0 ^ 2 ' c iBxp-\~— c d Gz P -\r~2 

(97) 


Here, the moment M Y of course applies to the case 
of the similar-section wing, as given bv equation 
(92). 


The lift of the cambered-span wing can also be 
expressed in coefficient form 

Ci-yl— (99) 

IpF-S 

and since C L is related to the section force-coeffi- 
cient distribution Ci(s) and chord distribution 
c(s) by 

Ci=£ _1 f / CiC cos rds (100) 

J -s t 

the variation of M P and hence C m ,p with C L can 
be determined. 

It should perhaps be emphasized that the 
function C m ,p(C L ) does not represent the variation 
of C m ,p as the angle of attack of a particular wing 
is varied, as with conventional airfoils. Rather, 
it gives the pitching-moment-coefficient variation 
with the optimum wing form corresponding to the 
maximum wing LJD at each value of C L . That is, 
the physical wing form is different for each value 
of C L . 

VARIATION OF C m WITH Cl FOR SIMILAR-PROFILE WING 

The pitching-moment coefficient of the wing 
constructed of similar section profiles all operating 
at the same value of c t is also given by equation 
(98), but with the value of M P being obtained 
from equation (97). In this special case, however, 
equation (100) reduces to the simple form 

r L =- (ioi) 

m 

where m is a constant (eq. (27)). Thus, the 
individual terms of equation (97) can be directly 
evaluated for each value of C L , using the corre- 
sponding section force coefficient Ci=mC L . 


VARIATION OF C m , P WITH Cl FOR GENERAL WING FORM 

Equation (90), in conjunction with equation 
(85), gives the wing pitching moment M P about 
any arbitrary axis. In order to express this 
equation in dimensionless form, a pitching- 


LOCUS OF TRIM AXIS 

As can be seen from equation (90) for the most 
general form of the pitching moment, there are 
an infinite number of axis locations for wdiich 
M P will be equal to zero. An axis for which 
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Afp= 0 will be defined as a trim axis. A knowledge 
of the locus of such axes is obviously important 
for determining aircraft component arrangements 
which will possess satisfactory trim properties at 
cruise flight conditions. 

For a given set of cruise flight requirements 
and the corresponding optimum wing form, the 
value of M y (eq. (90)) can be calculated and the 
values of the two integral factors determined. 
Then, the trim condition M P = 0 leads to the 
equation 


Z T = 


if ‘ 

J -s t 


M Y -\-x T q I CiC cos rds 


<l j ' cds 


( 102 ) 


which defines the locus of the trim axis, z T =z T (x T ) . 
Since M Y and the integrals are constants in this 
equation, the z T value for any x T location is thus 
determined for the trim condition. 

For the case of a wing constructed of similar 
profiles operating at the same value of c h the locus 
z t — z t (x t ) is obtained from equation (97), 


Zt — 



qc 0 CiB 



b I Wo T> b 
2 + V ClB 2 


r] 


(103) 


By use of t lie equations developed herein, the 
wing pitching moment of any optimally loaded 
oambered-span airfoil can be calculated for the 
design flight condition. In the general case, the 
optimum nondimensional circulation loading T/r o 
may have to be determined by 1 he electrical 
analog method and hence must be expressed in 
graphical form. ( Consequently, the ensuing de- 
sign analysis and moment determinations will 
also have to be carried out graphically. Even in 
the particular cases where F/T 0 can be obtained 
analytically by conformal transformations, the 
resulting functions may result in extremely 
complex integrals, so that machine solutions are 
necessary. 

In the particular case of the similar-profile 
wing, the calculation of the pitching-moment 
variation with lift coefficient is considerably 
simplified because most of the integrals are then 
constants. 

The locus equations (102) and (103) can be 
used to determine the proper location of the 


center of gravity of a complete aircraft configura- 
tion so as to optimize the longitudinal trim 
requirements at cruise. 

CONCLUDING REMARKS 

General relations needed for the design of 
arbitrary cambered-span airfoils which will possess 
the theoretical minimum induced drag for specified 
flight conditions have been developed. These 
relations, however, allow the determination of the 
optimum wing form not only for minimum in- 
duced drag, but also for the maximum attainable 
value of lift-drag ratio by optimizing the chord 
distribution of the wing. 

The procedures developed can be used for the 
determination of the optimum wing design when 
a specific spanwise camber line and section pro- 
file are given. The camber-line shape and wing 
section profile that will be best for a particular set 
of flight requirements depend, of course, upon the 
specific mission involved, and only by a series of 
comparative designs can the best overall wing 
form be determined in any particular case. The 
design procedures presented, however, are in a 
form which allows such efficiency comparisons to 
be easily made. 

It should perhaps be explicitly emphasized that 
cambered-span wings will not in all cases possess 
superior aerodynamic efficiency compared with 
the optimum flat wing of equal projected span. 
The relative efficiencies of flat and curved wings 
will depend critically upon the ratio of the induced 
drag reduction to the increased profile drag of 
the cambered wing, and t lie magnitude of tins 
ratio can be determined only by carrying out a 
series of comparative design analyses according 
to the methods presented herein. 

The design procedure presented does not specifi- 
cally include such possible effects as interference 
and tip separation drags. The method is based 
on the assumption that the two-dimensional sec- 
tional force coefficients are closely approximated 
by the profiles of the three-dimensional wings. 
For large span wings or wings with moderate 
spanwise camber, the procedure outlined should 
be quite valid within the limitations of linear 
lifting-line theory. The procedure also includes 
the assumption that the section profiles selected 
for wing construction have a sufficiently large 
minimum thickness that they will provide ample 
room for housing the wing spar structure. 
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The results of the illustrative comparison of a 
cambered wing with a flat wing of equal span 
(although neither wing has been optimized with 
respect to maximum lift-drag ratio) indicate that 
gains in operational efficiency can be secured with 
cambered-span wings, as compared with the effi- 
ciency of conventional wings currently in use, 
when use is made of special profiles to minimize 
the wing profile drag. The magnitudes of the 
gains attainable with cambered-span wings over 


optimum flat wings will be smaller at low lift 
coefficients due to the predominance of the profile 
drag. At higher lift coefficients, however, where 
the induced drag becomes a significant factor in 
wing efficiency, the higher effective aspect ratio 
of cambered wings becomes very important for 
wings of limited span. 

Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Station, Hampton, Va., June -5, 1962. 


APPENDIX A 


THE RELATIONSHIP OF LIFT TO INDUCED DRAG FOR OPTIMALLY LOADED AIRFOILS 


It can be shown (ref. 1) that for any optimally 
loaded lifting line of projected span b f and pro- 
ducing a lift force L , the induced drag is given by 




KF 


(Al) 


where q is the free-stream dynamic pressure and/ 
is a constant that depends only upon the shape of 
the spanwise camber line z(y). 

For a straight lifting-line segment, the circula- 
tion distribution is elliptical when the line 
optimally loaded and 


is 


irq 


CpJ 


Therefore 


/ (straight line) — -=0.318 


7 r 


(A2) 


(AS) 


For a cambered-span airfoil (or any other non- 
planar lifting system), the induced drag is 

a= (RV0<z[f] (A4) 


and 


/ (cambered (AS) 

From equation (4), when ^=1.0 (corresponding 
to the condition where the projected spans of the 
straight and cambered wings are equal), 

Jc = 7r ~ 1 BN A (A6) 

and the / value for any cambered form, when 
^=1.0 is used to determine k, is then given by 

1 0.318 


/ (cambered line): 


7T k 


(A7) 


When the cambered span will possess less 

induced drag than the optimally loaded flat span. 
Thus, for wings of equal projected span, the 
induced drag depends only upon the total lift L 
and spatial distribution of the wake vorticity 
(whose effect is measured by/), and is independent 
of wing area S or the associated aspect ratio 
/ b' 2 \ 

Al=^g ) so long as the wing loading is optimum. 

The need to introduce a reference area S and 
aspect ratio A arises in defining coefficients. 
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